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Abstract 

We perform a detailed analysis on the scaling properties of the total 
7*p cross section, <t 7 » p . We write the cross section as a product of two 
functions W and V representing, respectively, the dynamical degrees of 
freedom and the contribution from the valence partons. Analyzing data 
from HERA and fixed target experiments we find that V is nearly in- 
dependent of Q 2 and concentrated at large x, while W carries all the 
information on the Q 2 evolution of 7*p. We define the reduced cross sec- 
tion o" 7 * p = W = <j~,*p/V, and show that it is very close to a generalized 
homogeneous function. This property gives rise to geometric scaling not 
only for small x, but for all the current measured kinematic plane. As a 
consequence of our Ansatz we also obtain a compact parameterization of 
<7 7 *p describing all data above Q 2 = 1 GeV 2 . 
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1 Introduction 



It has been found that for low values of the Bjorken variable x, x < 0.01, the 
total 7*p cross section, <7 7 * p (a;, Q 2 ), extracted from lepton-hadron scattering 
presents the property of geometric scaling ^ This property permits to 
write the cross section as a function of only one variable, r, called the scaling 
variable, which is the product of two functions, one depending only on Q 2 and 
the other only on x. It has been suggested that, for ct 7 » p , t is given by Q 2 /Q 2 
with Q 2 — Ql(x) known as the saturation scale. Geometric scaling has also been 
observed in eA reactions , inclusive charm production 4 and nucleus-nucleus 
collisions 

The observation that <t 7 * p grows quite rapidly at low x and that this behavior 
can not continue indefinitely without violating the unitarity of the cross section 
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led to the proposal of nonlinear QCD equations containing saturation |S1[7||S1|5J 
1101111) . One of the features of this type of equations is the introduction of a scale, 
Q 2 , to signal the onset of saturation effects. Much of the excitement and advance 
in the understanding of perturbative QCD at low x in recent years comes from 
the discovery that some saturation equations imply geometric scaling at the 
saturation scale [T21 IT31 HI ITS] . 

Much recent theoretical work has been devoted to find both, the region 
where geometric scaling is valid and the functional form of the violations to 
geometric scaling above the saturation scale. Studies based on the BFKL equa- 
tion supplemented with specific boundary conditions have found that there 
is a window of phase space above Q 2 and below a Q^ ax on which geometric 
scaling is valid. For current accessible energies, Q^ ax is of the order of 100 
GeV 2 |171 1181 H5| . More recently, there have been indications that, in a more 
general nonlinear equation, geometric scaling is strongly violated [201 - 

As the concept of geometric scaling has been linked to saturation, none of 
these investigations expect geometric scaling to be valid at large x where the 
density of partons is very small. 

Here we study in detail the scaling properties of the total 7*p cross section 
and find that geometric scaling is related to the fact that at small 
very close to a homogeneous function, specifically a power law in both, x and 
Q 2 . We show that it is possible to define a reduced cross section, called <t 7 * p 
in the following, which isolates this power law behavior not only for the small, 
but also for the large x region and thus shows geometric scaling in the complete 
kinematic plane. 

This document is organized as follows. In the next section, we show that it 
is possible to isolate the power law behavior in x of cr 7 * p for all values of Q 2 and 
define the reduced cross section <t 7 * p . In Section|3]we study the behavior of ct 7 * p 
and show that it is very close to a generalized homogeneous function. In Section 
0]we discuss the implications of our findings regarding saturation and geometric 
scaling. We also present a compact parametrization of cr 7 » p which describes 
all data above Q 2 = 1 GeV 2 . Finally, in Section we briefly summarize our 
findings and present our conclusions. 



2 Analysis of the x and Q 2 dependence of cr 7 * p 

First, we turn to the behavior of cr 7 * p at small x. It is know, that the experi- 
mental data at small x can be described at each value of Q 2 by a power law in 
x [22] • From the point of view of theory, this behavior is expected, for Q 2 big 
enough to justify the use of pQCD, from both, the DLLA approximation if 
the starting Q 2 value for the QCD evolution is taken sufficiently small (see for 
example [21]), and from the BFKL evolution. 

In fact, this behavior is also seen in studies of geometric scaling above the 
saturation scale. For concreteness we use the saturation scale Q 2 {x) ~ i _Agbw , 
with A GBW = 0.288, as defined by K. Golec-Biernat and M. Wiisthoff |3T]. 
Stasto et al. PQ found that, for r = Q 2 /Q 2 {x) > 1, cj 7 * p (t) ~ 1/r suggesting a 
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power law behavior of function of x for constant values of Q 2 . 

In summary, it is expected from theory and phenomenology, and confirmed 
by experiment that, at small x, er 7 » p behaves as a power law for fixed values of 
Q 2 . Following these results, we propose to write the total 7*p cross section at 
each Q 2 value as the product of two functions W and V: 

a rp = WV (1) 

where 

W = Nx~ x , (2) 
and we use x instead of x-Bjorken Both are related through 

(3) 

with rrif = 140 MeV. Note that x and x are substantially different only at very 
low Q 2 . Note also that in equation J2J) both, the normalization, N, and the 
exponent of the power law, A, may be different, for different values of Q 2 ; i.e., 
N = N(Q 2 ) and A = A(Q 2 ). 

As we are interested in isolating the scaling behavior of the cross section we 
require from V to be approximately constant at small x, so that it does not 
alter the physics embodied in W, and that it describes the cross section at large 
x. It turns out that these requirements are very closed to those expected from 
a valence distribution, so we tried the following functional form inspired in the 
work of [251 : 

/ {x~x ) 2 \ f l-(x-x ) \ 
= CX P { £5- J Crf { ^ J ' (4) 

where the Gaussian represents the valence distribution of the proton and the 
error function takes care of enforcing the different kinematic constrains. 

We proceed now to test the Ansatz embodied in equations (P^J . We use data 
from fixed target [27j and HERA [2111221 experiments. We chose Q 2 values such 
that there are measurements from both, fixed target and HERA experiments, 
either at the same or at very similar Q 2 . Normally, these Q 2 values are slightly 
different. We correct the fixed target data to the HERA values using the HI 
PDF 2000 fit [2H1 - In most cases the correction factors are at the per mil level 
with a few cases at the one and two percent level. We then fit the data to 
the functional form of equation for each Q 2 value. As an example, Figure 
n shows data from HERA and fixed target experiments at Q 2 =12 GeV 2 and 
Q 2 =120 GeV 2 along with the result of the fit. 

The data in Figure ^ is very well described by equation QJ. Furthermore, 
the same can be said at each value of Q 2 where there are measured points from 
both, fixed target and HERA experiments. Remarkably, we find that xq and a 2 
from equation (@J do not depend on Q 2 . They take the values xo = 0.27 and 
cr 2 = 0.036. 
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Figure 1: The total 7*p cross section is shown as a function of x for Q 2 =12 
GeV 2 and Q 2 =120 GeV 2 . The points are experimental data from fixed target 
and HERA experiments and the line is a fit to the form of equations 
A power law is clearly seen at low values of x, while the structure at high x 
corresponds to a Gaussian-like distribution. 
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The fact that V is independent of Q 2 implies that all the QCD evolution 
in Q 2 of the cross section is entirely contained in W. Thus, we define, in the 
complete kinematic plane, the reduced cross section <7 7 . p as 

<V P = W(x,Q 2 ) =a rp /V(x). (5) 



3 Analysis of the behavior of <t 7 * p 

We turn now to study the reduced cross section. Specifically, we study the 
dependence on Q 2 of both N and A. We use all HERA data [211123 to fit <7 7 . p 
for fixed values of Q 2 . We do not use fixed target data at this stage of the 
analysis, because those data points are concentrated at high x and thus, they 
do not have a lever arm long enough to determine accurately the power law 
parameters. Their influence has already been taken into account through V(x). 

We use 40 different experimental Q 2 values, ranging from 0.15 to 8000 GeV 2 , 
with enough data points in x to perform the fit. The average number of points 
for each fit was 8, ranging from 5 to 12. At each value of Q 2 , equation |J2J 
provides an excellent description of data. The results for N and A for each 
individual fit are quite precise and provide a clear picture of their Q 2 dependence 
as shown in Figure |21 

We observe a dramatic change in the behavior of both functions, N(Q 2 ) 
and X(Q 2 ), when the virtuality of the photon approaches from above the region 
below 1 GeV 2 . Here A is, as expected, very similar to the one found with the 
Donnachie-Landshoff parametrization |3U) : while the normalization appears to 
saturate. Furthermore, the value of N and A for the Q 2 data below 1 GeV 2 are 
almost constant in comparison to the steep dependence of these functions above 
1 GeV 2 . 

Above Q 2 w 1 GeV 2 , A(Q 2 ) can be described with the following functional 
form 

A(Q 2 ) =alog 10 (Q 2 /A 2 ), (6) 
while N(Q 2 ) behaves as the power law 

/ o a\-(i+«0 

"(O- ■ CD 

where Qq is taken as 1 GeV 2 . 

A fit in the intermediate Q 2 region to the data plotted in Figure |21 yields 
= 41.1 ± 1.5 fib, e = 0.104 ± 0.007, a = 0.135 ± 0.003 and A 2 = 0.17 ± 0.03 
GeV 2 . The quality of the fits is x 2 /dof = 0.66 and yf/dof = 0.33 for N{Q 2 ) 
and A(Q 2 ) respectively. Note that the points at the largest Q 2 were not taken 
into account, because they have big fluctuations due to the limited statistics of 
data. 
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Figure 2: The Q 2 dependence of the normalization iV and the exponent A of 
<r 7 »p = Nx~ x , extracted from fits of HERA data to the reduced cross section 
tr 7 » p at fixed values of Q 2 . The solid lines are fits to equations © and (JJJ) in 
the intermediate Q 2 range given by the empty bullets in the figure. 
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In summary, for Q 2 values below 1 GeV 2 the Q 2 dependence of W is very 
soft and the function depends mainly on one variable: x. Above 1 GeV 2 W is 
very well described by the following functional form: 

\QoJ 

4 Discussion 

Scaling. The Q 2 dependence of A is only logarithmic. Consider first the case 
where the exponent of x is a constant A — Ao- Naming Wo the function W with 
Ao, equation (JHJ) would be of the form 

W (x,Q 2 ) =k(Q 2 )-^x- x °, (9) 

where k is just the normalization and where we have reverted to x instead of x 
for simplicity, because both are equal in this kinematic domain. 

In this case Wq is a generalized homogeneous function, which, as can easily 
be demonstrated, implies that for all t real and bigger than zero the following 
equation is valid: 

W (t- 1/X °x, t 1/{1+e) Q 2 ) = W (x, Q 2 ). (10) 
In particular, it is also valid for t — x x °: 

W Q \ t=x > B (l,x x °^ 1+ ^Q 2 ) = W (t ) = W (x,Q 2 ); (11) 

i.e., Wo exhibits exact scaling with the scaling variable given by to — x x "^ 1+ ^Q 2 . 

Now, we turn to the real case where A depends on Q 2 . Using the same 
t = x x ° one obtains 

W(r ) = W(x,Q 2 )x- 5 , (12) 

where 

5 = A -alog 10 (Q 2 /A 2 ). (13) 

In this case the scaling is broken when S, the exponent of x in the RHS of 
equation Ijl2(l is different from zero. With the help of Figure note that for 
values of Q 2 between 1 GeV 2 and the upper limit of Q 2 used in pp , Q 2 = 450 
GeV 2 (implicit for x < 0.01), the average value of A(Q 2 ) is close to the value 
of Aqbw- Therefore, for A = Aqbw the exponent of x in the RHS of equation 
<|12[) takes its smallest values. Furthermore, one can check that for any given 
Q 2 in 1 < Q 2 < 450 GeV 2 , the range covered in x (x < 0.01) is limited in 
the current measured kinematic plane to an average of one order of magnitude. 
That means that for the given Q 2 the departure from scaling for all x points 
is similar and numerically small. As in the small x region er 7 » p ~ <r 7 * p , this 
explains the approximated geometric scaling observed in pQ. 



7 




T T 

Figure 3: The left panel shows the reduced 7*p cross section, er 7 * p , as a function 
of the scaling variable r = Q 2 /Q 2 S [H] for all data in [27J US EH The geometric 
scaling behavior is clearly seen for all values of r. The right panel shows the 
scaling combination of equation (|12(l . Data points in the small x region, x < 0.01 
are shown as full bullets to have a comparison with previous studies £Q. 

This argument is valid in the complete kinematic plane because, as pointed 
out before, for Q 2 below 1 GeV 2 , W is already approximately a function of 
only one variable, namely x, and the data shows a smooth matching of the two 
behaviors. Thus, CT 7 *p should display approximated geometric scaling for all 
values of x and Q 2 . 

The left panel of Figure [3] shows that this is indeed the case, while the right 
panel shows the scaling combination of equation 1)12(1 . Note that this latter case 
can not be considered properly as scaling and that it is shown for illustration 
purposes only. 

The data span some six orders of magnitude in x, and another six in Q 2 . 
With the addition of data with x > 0.01 the geometric scaling behavior is ex- 
tended two orders of magnitude in r. It is quite interesting to compare Figure[3] 
with Figure 01 The latter figure contains all data points above Q 2 = 1 GeV 2 
before the data collapse produced by the transformation to the scaling variable 
r. The comparison of both figures shows that the collapse of all data in a single 
line is not a trivial fact. 

Power laws, scaling and critical phenomena. It has to be emphasized 
that the origin of scaling within this approach is the fact that W is very close to 
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a generalized homogeneous function. This fact is valid even at very large values 
of x, where one would not necessarily expect saturation effects to be present. 
But it does not exclude the possibility that the mechanism which gives rise to 
the power law behavior of W is also linked to saturation. 

In this context it is interesting to note that scaling and its relation to power 
laws has been widely discussed in relation to critical phenomena. In particular 
it has been found that under some conditions the presence of a renormalization 
group equation helps to explain the appearance of power laws, of its associated 
scaling and permits also to explain and numerically estimate the appearance of 
scaling violations (see for example 152] and references therein) . 

It is interesting to note that the power law behavior of the total 7*p cross 
section is generated by the branching process embodied in QCD evolution equa- 
tions which are in fact a type of renormalization group equations. Also the case 
of saturation has been cast, within the Color Glass Condensate approach in 
the form of renormalization group equations. It is clear then, that the subject 
of finding a deeper understanding of the relation between renormalization group 
equations and the emergence of power laws in pQCD deserves further studies. 

A parametrization of er 7 . p above Q 2 = 1 GeV 2 . Note that as a conse- 
quence of the description of cr 7 * p given by equation (JIJ we also have a simple 
six parameter description of the total 7*p cross section for all Q 2 values above 
1 GeV 2 : 



/ n 2\ -(i+O 

W ° y (14) 
exp ( _ £Z£2T\ erf^- ( "° r 



4ct 2 J V 2(7 

Equation l|14fl is compared to data in Figure0]using the parameters obtained 
from the fit to Figure|3| For these parameters the x 2 jdof obtained for Q 2 > 1 
GeV 2 is x 2 /dof = 0.77. 



5 Summary and conclusions 

We have shown that ct 7 * p can be separated as a product of a power law func- 
tion, W, carrying all the information on its Q 2 evolution and a Gaussian-like 
distribution, V, which depends only on x. As a consequence of this Ansatz we 
obtain a six parameter description of all cr 7 » p data above Q 2 ~ \ GeV 2 , where 
each parameter has a natural physical interpretation. 

We define a reduced cross section, W = ct 7 * p , and find that it is very close 
to a generalized homogeneous function over all the measured kinematic plane. 
This property is found to be responsible for the geometric scaling behavior of 
er 7 » p in the small x region and of <t 7 » p in the complete kinematic plane. These 
results show that the emergence of geometric scaling is not necessary related to 
saturation and open up interesting possibilities for further studies of the relation 
between evolution equations and the appearance of scaling behavior in QCD. 
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Figure 4: The total 7*p cross section, ct 7 * p , is shown as a function of x for 
different fixed values of Q 2 going from Q 2 = 1.2 GeV 2 to Q 2 = 20000 GeV 2 . 
The bullets are the experimental data points from HERA, while the triangles 
are from fixed target experiments. The alternation of full and empty symbols 
is just to get a clear display of data. The lines are the result of equation (|14|) . 
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